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A quantum Langevin equation for the amplitudes of the collective modes in Bose-Einstein con- 
densate is derived. The collective modes are coupled to a thermal reservoir of quasi-particles, whose 
, elimination leads to the quantum Langevin equation. The dissipation rates are determined via the 

correlation function of the fluctuating force and are evaluated in the local-density approximation 
for the spectrum of quasi-particles and the Thomas-Fermi approximation for the condensate. I take 
great pleasure in dedicating this paper to Gregoire Nicolis on the occasion of his sixtieth birthday. 
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I. INTRODUCTION 

o 

lyj ' The realization of Bose-Einstein condensates of very rarefied evaporatively cooled gases of alkali atoms in magnetic 

, traps 1^ offers the unique possibility to test ab initio many-body theories in the laboratory . One very fertile field 

' has been the experimental and theoretical investigation of collective modes of the condensates, both for zero and 

^ finite temperatures. For reviews of the experimental and theoretical work see and |^,^ respectively. As opposed to 

\^ conventional superfluids like He-II j^J^ in the new systems the collision-less regime is very naturally realized. In this 
regime the dominant damping mechanism for collective modes is Landau-damping, whose temperature dependence 

Q I in spatially homogeneous condensates in the regime ksT large compared to the chemical potential ^ has first been 

Q studied by Szepfalusy and Kondor Recent investigations , though more exact, led to similar 
results, differing by a prefactor close to 1 for the damping rate. For condensates in traps Landau damping of low-lying 

^ ' modes is more difficult to calculate, and additional approximations are needed to cope with the fact that momentum 

\ is not conserved in a trap. The damping rate of collective modes in traps has been calculated in Q using the local 

• density approximation and, in addition, a classical approximation for the correlation function whose Fourier-transform 

' determines the cross-section of Landau-scattering. For the isotropic breathing mode in isotropic traps the Landau- 

, damping has been calculated numerically by evaluating the coupling to a great number of discrete quasi-particle 

■ modes and subsequently introducing some smoothing. The quasi-continuum coupled to the collective mode under 
, study was displayed explicitly in this work. Theories using an extension of the approach of |^ via the dielectric 
' formalism and an approach via a time-dependent mean field scheme ||T7| have also been given. 

— ] In the present paper a very direct approach |18| to the dissipative equilibrium and non-equilibrium dynamics of 

■ collective modes in trapped Bose-Einstein condensates via quantum Langevin equations is put forward. Because of 
' the discreteness of the mode-spectrum in traps the problem is formally similar to the quantum-optical problem of 

I discrete modes in a laser, for which the formulation in terms of quantum Langevin equations has been very useful 

Ph , In the next section the microscopic description of a trapped Bosc-Einstcin condensed gas is briefly set up. Then 
^ ' we recall the basics of the quantum Langevin equation of a boson mode. The derivation of the quantum Langevin 
equation of a collective mode follows. The damping rates are then evaluated in the local density and the Thomas-Fermi 
approximation. The last section contains a discussion of our results. 



C3 ■ II. MICROSCOPIC EQUATIONS OF MOTION 

The weakly interacting Bose-gas in a trap in standard notation is described by the Hamiltonian 

H = fd^x^P+l - —S/^ + V{x)- (1) 
J I 2m 2 J 

The presence of a Bose-Einstein condensate means that many (Nq 3> 1) particles occupy a single normalized mode of a 
macroscopic classical matter wave, determined as the mode of lowest energy of the classical Hamiltonian corresponding 
to eq. ^ . It satisfies a classical wave equation, the Gross-Pitaevskii equation , which we take in an extension defined 



1 



by the so-called Popov-approximation |2|] , including the interaction of the condensate with the density n' of thermal 
atoms, but neglecting its interaction with the pair amplitude ~ V'o of thermal particles 

- (;iV2m)VVo + {V{x) + Ua{Na\Mx)? + 2n'(a;))) = mV-o- (2) 

For given TVq the chemical potential follows by imposing the normalization condition on ipQ. 

The presence of the highly occupied condensate mode makes the decomposition of the Heisenberg field-operator 

V'(a;, t) — (^y/NoijjQ{x) + ijj'{x, t)) exp {—i^t/Ti) exp (3) 

useful. \/Nq exp is the complex amplitude of the classical condensate mode in equilibrium. tp'{x,t) is the field 
operator for the particles outside the condensate. We shall neglect fluctuations of the number of atoms in the 
condensate and also fluctuations of the phase of the condensate, which can be shown to occur on a time-scale much 
longer than the relaxation-time of the collective modes p^ . 

The Hamiltonian splits up according to H = Hq + Hi + H2 + -ffs + H^, with a c-number term Hq which need not 
concern us here, and 



Uo{No\ijJo\'^ + 2n'))?/;Q'0' -I- {hermitian conjugate) ) 
Hz = Uo\/l% / d?x (^tpQ{'4>''^i/j' — 2n')'0' -t- {hermitian conjugate)^ 



H^^^ J d^xli^'+i^'+i^'i}' - 4n'iP'+iP'). 

The splitting is here done in such a way that the term Hi vanishes due to eq.(|^) and the part H2 describes the 
linearized quantum excitations around the solution of (|^) . H2 is diagonalized by introducing quasi-particle operators 
d,y, a'^ by the standard Bogoliubov transformation 

Tp'ix^t) = ^ {uy{x)au{t) + vl{x)cil {t)) , 

where Ui,^ Vi, satisfy the usual Bogohubov-De Gennes equations 

( -i;:^^ + U,s{x)-hu, K{x) \ (u,{x) 



2m 

K*{x) -|iv2 -t- U,s{x) + hco, ) \vu{x) 



0, (4) 



with the abbreviations 



U,s{x) = V{x) + 2t/o(iVo|Vo(a;)P + n'{x)) - 

K{x) = NoUoMxf ■ (5) 

Eq. is consistent with the ortho-normality conditions / d^x{uyU*^ — = and / d'^r{u*Vfj, — = 0, 

which guarantee the Bose commutation relations of the a^, a^. The decomposition of ^ and tp' together imply that 
iV = iVo + / n'ix)d^x with n' = M^^l^ + I^mP) + 

Within Bogoliubov (-Popov) theory the terms H3, H4 of the total Hamiltonian are neglected and the quasi-particle 
operators 6;^ in the Heisenberg-picture obey the Heisenberg equations of motion dj/ = —iuj^a^. In this approximation 
the collective modes and the quasi-particles have infinite lifetime. In reality, however, the lifetime will be limited by 
the scattering of quasi-particles in any given mode v with other quasi-particles from the thermal reservoir, which is 
described by H^ and H4. One way to describe this is the quantum Langevin equation. 
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III. QUANTUM LANGEVIN-EQUATION OF A HARMONIC OSCILLATOR 



Let us recall here briefly the quantum Langevin equation, in Markoff approximation, of a harmonic oscillator as it 
is commonly used in quantum optics [^|o). For a detailed discussion of its derivation I refer to |2^. Eq. (3.4.63) of 
that reference states the quantum Langevin equation in resonance or 'rotating wave' approximation for a harmonic 
oscillator, described by the Bose operators 6,0+, in interaction with a thermal reservoir at temperature T. It takes 
the form 

a(t) = -ina{t) - 7a(i) + ^(t) (6) 

where fl is the frequency of the oscillator including a frequency shift due to the oscillator's coupling to a heat reservoir, 
7 is the damping rate, and ^{t) is a Gaussian noise-operator. In Markoff approximation it has the correlation functions 

exp(Ml/ kb^ ) — 1 

ensuring the correct normally ordered expectation values in equilibrium, and 

(m,i+{t')]) = 2^5{t-t') (8) 

ensuring the correct Bose commutation relations of a(i),a+(t) for all times. The fluctuation-dissipation relation 
therefore permits us to infer the properties of ^(t) if the coefficient of the dissipative term is known. Alternatively we 
can infer the dissipation rate 7 from a microscopic expression for ^(t) either by using (|^) or, alternatively, (^). 



IV. QUANTUM LANGEVIN-EQUATION OF COLLECTIVE MODES 

We shall here confine our attention to the dynamics of the low-lying collective modes in the collision-less regime. 
The interaction of the collective modes with the thermal quasi-particles is described by the Hamiltonian + not 
yet contained in the Bogoliubov(-Popov) approximation. Because it contains the large factor \/ Nq the contribution 
H3 dominates over H4 and the latter can be neglected in the following. Inserting the Bogoliubov transformation 

in H3 and going to the interaction picture with respect to the unperturbed Bogoliubov-Popov Hamiltonian, H3 in 
interaction representation takes the form 

^3 = ^ E + J*)a+d.d^ exp[zK - - ^^)t] 

+ {hermitian conjugate)} + [nonresonant terms). (9) 

where = l/(eKp{huj^/kBT) — 1) is the thermal number of quasi-particles at frequency uj^. 

Nonresonant terms, in which the frequencies of the quasi-particles cannot add up to zero, have not been written 
out explicitly, because later-on we shall restrict ourselves to the resonance or rotating wave approximation in which 
they don't contribute. The relevant matrix elements M'^^^M^^-' are 

Mill, - 2Uo J d^xiJou:{v;v, + ^u;u,) + fi). (10) 

Mk]1^ describes a scattering process in which one atom is scattered out of the condensate by the absorption of the 

two quasi-particles out of and the emission of the new quasi-particle k into the thermal bath. Likewise Mlfl^^. 
describes a scattering process where an incoming thermalquasiparticle k is absorbed, again an atom is kicked out of 
the condensate, and two quasi-particles fi are emitted into the thermal bath. The scattering amplitudes for both 
processes are linearly superposed due to the phase-coherence of the condensate on the time-scale of the relaxation 
process induced by the scattering process. 

Taking into account the equations of motion ofce^{t) in the interaction picture a,y{t) = a,y{t) exp{—iuj^t) become 



3 



+2 [iM(%r + Mi^J,,] a^tttK exp[zK + - ^k)*]} ■ (11) 



If the back-action of the collective mode on the quasi-particle operators d^, in ( |llD can be ignored, the new term 
in this equation of motion acts like an effective random force operator. In the resonance approximation the average of 
this force operator vanishes. In addition it is white noise, in good approximation, if the frequencies uj,^ — lu^ — and 
lUk+uj^ — ijjy it contains form a closely spaced quasi-continuum near in a neighborhood which is broad compared to 
the resulting damping rate 7^. For an explicit display of this quasi-continuum in a concrete example see W%. In as 
much as this condition is satisfied for large condensates the Markoff-assumption made earlier is justified. All terms 
in the fluctuating force term not containing frequencies near frequency are non-resonant and can be omitted in 
comparison with resonant terms. 

As we recalled in the previous section, the noise term is always accompanied by a dissipative term, and, due to 
the Kramers-Kronig relation, also by a frequency shift. Thus the complete quantum Langevin equation in resonance 
approximation and in Markoff approximation must take the form of (^) 

OLy ^ -iiyiy + &^)a^ - ^^a^ + iy(t) (12) 



where is given by the second term in ( Jl] 

The damping rates will be derived below, but we can also simply use and represent them in the form 

7.= dt{\Ut)Xtm)- (13) 

Evaluating the commutator, taking the thermal expectation value, and performing the time integral in (|l^) we obtain 

+ + {M^'Jjrin. + ^)<5K +^,^ ^.)} (14) 

The first term describes Landau-damping of the mode v by scattering a quasi-particle from mode fi to mode k and is 
equivalent to a result derived in jl^] by the golden rule. The second term in eq.(|lj) describes Beliaev damping, where 
the mode i' decays into two modes k, fi. It survives even for T ^ where — > for all modes. However, for low-lying 
modes in traps there are only very few modes, or no modes at all, into which decay under energy conservation can 
occur, and this contribution to the damping is then negligible. 

Let us now derive the dissipative term of the quantum Langevin equation. To this end we consider the equations 
of motion for ^(o^^a^) and ^(a^d^)! keeping again only the resonant terms. Integrating these equations over time 
from — 00 to t and inserting the result back into the equation of motion for ai, we obtain 




^° E + {Mi%r]&.{~^)M-^)exp[t{u;, - uj^-uj,)t] (15) 

KfJ, 

+2[{Mi%r + M(2) Jd+(-oo)d«(-oo) exp[*(c., +u;^- w„)t]} 

where the limit e — > +0 is implied. The second term on the right hand side is the fluctuating force term again, 
now more rigorously expressed in terms of the reservoir operators at the initial time at —00. Taking the limit with 
(e — iuj)~^ — !■ TrS{LL>) + iP/u>, where P/lo denotes the principal part under a frequency integral, we obtain the result 
( p^ ) for the damping rate and also the frequency shifts S^, in the quantum Langevin equation. They are given by the 
Kramers-Kronig relation 

= -ip / du.^ (16) 

where we defined ^{ijJv) — 7i/- 
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V. DAMPING RATES 



In the following we shall neglect the second term in (14), because as discussed it cannot contribute for low lying 
modes. Our goal in this section is the evaluation of the first term in ( [1^ ) in a well defined approximation, the 
local density and the Thomas-Fermi approximation. The local density approximation amounts to the treatment of 
the quasi-continuum of the spectrum of frequencies 0;^ — tj^ — uji, as a continuum whose density is given by the 



semiclassical mode-densities of the frequencies 



Why these frequencies lie much higher than the collective 



mode frequency uji, will become clear below. The Thomas-Fermi approximation applies to large condensates |2J] 
and amounts to neglecting the kinetic energy term in the Gross-Pitaevskii equation. The collective modes satisfy 
El, ^ fi — ?7o|V'o(0)P and can be represented as [Ej 



Vt,{x) = - 




(17) 



with JcPx\xiyix)\'^ = 1. The functions of the low-lying states Xi^ix) are known p5| , p6 28|2^,|2^ in the hydrodynamic 
(long-wavelength) and Thomas-Fermi approximation (and neglecting the influence of the thermal cloud which sits 
mainly outside the condensate and therefore has little influence on its collective excitations). In spatially isotropic 
parabolic traps they have the form |^5|| 



xAx) = -y^Pt^\x/rTF){x/rTF)iY,^ 



' TF 



The normalized polynomials P^^"'\x) are known explicitly p6 | 
The high- lying quasiparticle modes can be represented as ||4] 

„2 



with the local energies in Thomas-Fermi approximation 



,{e,^)Q{l-x/rTF) 



(18) 



^2E^pl/m 



(19) 



E,=E{p^,x) = W( 



2m 



\Uono{x)\r~U§nl{x)e{fi-V{x)) 



(20) 



and noix) = No\Mx)\^ = Wf/o)(l - E^i^^/4F?): 

and r!^p = y^2fi/ujf are the three main Thomas-Fermi radii. 
Let us consider now the Landau-damping of a low- lying phonon mode uj^. If the modes fi, k involved in the scattering 
process were also low-lying we could use eq. (|l7| ) and would obtain, with Ef^ = E^ + E"^ 



r(l) ^* 1 Ji r(2) "^^0 f ,3 I * * \ Ep^E^Efi 



(21) 



However, in the limit of low- lying modes where Ei,,E^,E^ <C UoUq this matrix element becomes very small, i.e. 
low-lying modes cannot significantly contribute to Landau damping of other low-lying modes. Therefore the relevant 
modes /x, k are in fact not low lying, local density approximation is applicable, and we can use eq.(p^ for their 
representation. The matrix-element for E^ << Efj_,Ei^ can be expanded in E^/p, to lowest order around E^^ = E^ 
and becomes then 



with 



J d^xxl{x)exp{i{p^~Pf,)-x))F{Ef„p^) 



FiE,,p,) 



pI iEl + {pl/2mf 
2mEM + {pl/2mf 



(22) 



(23) 
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It will be very convenient later to express the product xt{x)Xi^{^') by the associated Wigner- function via 

j^e-^''<'^-^^W.{-{x + x'),k) (24) 

In the following we denote 

{x + x')/2^ X x-x'-*r (25) 
The rate for Landau-damping can then be written as 

with 



The sums over the states /x and k are only symbolic, because in the local-density approximation the discrete states 
have been replaced by a continuum which is normalized on the (5-function. Concretely, under the integral over x the 
sums over the energy levels of the scattering quasi-particles /i, k are in the local density approximation replaced locally 
by classical phase-space averages for fixed E'^ , and final integration over E^ and E^ which takes automatically 
care of the normalization on the (5-function. Thus 

i:,^'. - C - /"^J -E(p„..)y^. (28) 

In the following ^nd ^^^^ interpreted according to eg. (p8|) . The spatial integration over r can be done 
and produces a momentum-conservation factor (2T:h^^ S^^^^ (p^ ^ ^ ^^)- Next the integrations over and iS^ 
contained in can be performed, which just cancel the (5-functions of overall momentum and energy conservation 
and replace everywhere else E^ ^ E^^ + Ey and ^ Pfj. + hk. Then the expression for 7^ is reduced to 



7, = CEl / d^x 



_fk_ 
(27r)3 



iy,(a;,fc) J dE^ 



1 4^6{E,-E{p,,x))^^&^5{E, + E,-E{\p, + hklx)) (29) 



Next the integration over the directions of relative to k is carried out by using up the second of the two 
i5-functions explicitly displayed in eq. (^9|). This produces a factor 2ii for the azimuthal angle, and a factor 
\dE{pfj^, x) / dp'^^\^^{2p^hk)~^ from the integration over cos 6* between -1 and 1, where is the angle between fc and p^. 
Finally the integration over the absolute value p^ is done using up the last (5-function, which picks out the £c-dependent 



momentum- value — \/2m^ ,Je^ + UQnQ{x) — UQno{x). leaving us with the expression 

CEl /• ,3 f d^k W,ix,k) f^^ F\E^,p^;^^) 



4^2^37 7(271)3 J '^4/9£(p1^^n2 • , 2 / ^ ^ ' 



2feBT 



We now have to face the difficulty to evaluate the conditional average of {Tik) ^ . The rigorous way to do this, which 
unfortunately leads to multiple integrals which are tedious to evaluate, is to invert (^) which yields 

d^fc W,{x,k) _ /• 3 xl{x + r/2)xAx-r/2) 



(27r)3 hk J 27r2?ir2 

A much simpler way consists in expressing the desired average by the local sound- velocity ■\//i/mc^(a;) defined by 
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(27r)3 fik \ m 



Xu{x)\'. (32) 



and estimating the dimensionless sound- velocity c^{x) semi-classically as c^{x) w -^/l — [x/rTpY with the geometrical 

mean Thomas- Fermi radius ttf = Of course the use of the semi-classical approximation for the low 

lying collective mode is highly questionable and cannot be quantitatively accurate. Still we may like to use it as a 
rough estimate in a case where an accurate evaluation is not required or too time consuming. Below we shall check 
this approximation in two cases, where it cannot be expected to be particularly good. 

We now introduce scaled variables = Ef^/{NoUoipQ{x)) and x = x/ttf, x' ~ x' /ttf with dimensionless 

mode- functions Xi'(^) = ^TEXvix)- Altogether, using (|3l|), we are left with the result 

with the dimensionless function 

. fd^, /.3,,X:(^)x.(^;)(l-(^ + ^0V4) 
J J v*^ ) 

■I I dE, I . . ^ I . I (34) 



{El + l)smh{^E,{l-{x + ir/i)) 
For z >> 1 the functions become linear in z = kgT/fi and reduce to 



J J {x-xy{i- {x + xy/4) 

The result for the spatially homogeneous case Q can be recovered from eq.(^ for fc^T >> ^ by using the scaled 
homogeneous condensate density 1 — — > 1, the phonon energy E^ = ^ and normalized plane waves to 
evaluate / d^x J d^x'xl{x)x^{x)/{x-xf = (27r2/A;^)^TO<I>2/2/x which, together with (HjH) yields 7^. = ^ak^,!^. 

For isotropic traps the asymptotic result ( |35| ) becomes 



^ ^ d(cos 6')d(cos 6'')d<?!'/'™" (cos 6')P™- (cos 61') exp(-iTO^0) 

(2i£' (cos 61 cos 61' -I- sine sin 6*' cos 0) - 2)2 - {i^ + i'"^ - 2)2 ^^^"^ 
-1 -1 

where the functions PP (cos6) are the associated Legendre functions appearing in the spherical harmonics. 
If instead of (31) we use (|3|) to evaluate the conditional average of (hk)^^ we obtain in place of ( |3^ ) 

Hna.^A^) - / di^=f=iPt^\,),^^r (37) 

E^ J V(l-a;2) 

As was already emphasized, this result can only serve as a rough estimate for (^). 

In the simplest case — 1,£^ — 0, which is the isotropic fundamental breathing mode, we have PQ(cos(e)) = 
l,Pi°Hx) = |\/7(l-fa;2). In this case ( pq ) can be reduced to the numerical evaluation of a two-dimensional integral 
and we obtain 

70,0 - 26.42..c.o(a'no(0))i/2M^. (38) 

A* 

Eq.(p7|) yields via elementary integration 70,0 ^ 27 .27 . .ujo{a'^no(0)y ksT / fi which agrees surprisingly well with the 
more rigorous result (^. Can this be considered typical? The answer is negative: 



6^^ ^^^^V!^^::~^^' / rf^^^i-f"-^(^)^^- 1 dx'x''pt"\i'y' 
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The simple result (37) lends itself to further evaluation for modes with 7^ 0. For the surface modes with 



= 0, 7^ we obtain the estimate 

,o,c - -oia^MO)r-^"-Vi.^-f^ (39) 

In this case a numerical comparison with the more accurate expression ( ^ ) for the case £^ = 2 shows that the latter 
is about 30 percent smaller, probably giving us a realistic impression of the accuracy of the approximation for Ciy{x). 
For larger values of and the accuracy of this estimate can be expected to improve. 



VI. DISCUSSION AND CONCLUSION 



In the present paper the many-body problem of collective modes in Bose-Einstein condensates in interaction with 
thermal quasi-particles was addressed by a method based on the equations of motion of the quasi-particle operators. 
This method leads directly to a quantum Langevin equation for the creation and annihilation operators of the collective 
modes, containing fluctuating force terms, a dissipation term, and a frequency shift term. These quantities are related 
by the fluctuation-dissipation relation and the Kramers-Kronig relation. Each part of the interaction-Hamiltonian 
beyond the unperturbed Bogoliubov-Popov Hamiltonian in principle gives rise to separate contributions to all three 
types of terms. We have here considered only the most important of these, namely the part of the interaction 
Hamiltonian giving rise to Landau-damping. 

Dissipation can arise only from energy conserving real processes, which is manifest by the appearance of the 
energy conserving (5-functions in the expressions for the damping rates. This means that only resonant processes can 
contribute to these rates. In finite systems like the trapped condensates this causes a problem, because there the 
mode spectrum is discrete, the spectrum of frequency differences lu^ — uj^^ ~ ujjj near is only a quasi-continuum, and 
the dissipation rates in a strict sense have to vanish. In other words, in a strict sense, what is seen as dissipation 
can only be a 'short-time' effect; waiting for a sufhciently long time interval on the order of the inverse spacing of 
the quasi-continuum, revivals would have to appear. These will not be seen, however, at least in large condensates 
to which the local density and Thomas-Fermi approximation can be applied, because not only the energy stored 
in the collective mode but also the thermal energy of the system is available to bring into play a large number of 
modes which will lead to an irretrievable dissipation of the energy over many degrees of freedom. Therefore it is 
reasonable in such cases, if not required, to eliminate all recurrence effects, replacing the quasi-continuum by a true 
continuum, which is what the local density approximation does. Using this device we have arrived at definite results 
for the temperature-dependent damping rates of any collective mode, in an isotropic trap, which can be evaluated by 
computing numerically a multi-dimensional definite integral, e.g. by a Monte-Carlo routine. 

The different pieces of the perturbation Hamiltonian also each give rise to frequency shifts. These are generated 
by virtual processes which do not require energy conservation, i.e. resonance. However the effect of the non-resonant 
processes is suppressed by corresponding energy-denominators and small. Here we have limited our considerations 
only to those processes which can also become resonant. We have here not evaluated the frequency shifts further 
using the local density approximation as we have done for the damping rates. 

Experimental results for temperature dependent damping rates and frequency shifts have been obtained for 
anisotropic traps only |^], and we therefore refrain from a cornparison with our explicit results for isotropic traps. 
Detailed comparisons have been made for anisotropic traps in [p^ , where heavier and more powerful formalisms 
were brought to bear together with a stronger reliance on numerical work. 

The goal here has been more modest, namely to use a minimum amount of numerical work and to apply the direct 
and intuitive quantum Langevin approach to the fluctuations, damping rates and frequency shifts of collective modes 
in spatially inhomogeneous trapped Bose-Einstein condensates. 
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